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1. Introduction 

For quark confinement, the most well known and natural criterion is the Wilson criterion: 
the area law for the Wilson loop average, which means a linear static potential between a pair of 
quark and antiquark. The Wilson criterion for quark confinement is a gauge-invariant statement, 
which is independent from the gauge fixing chosen in quantizing the Yang-Mills theory. For gluon 
confinement, such criterion is not yet known as far as I know. 

For color confinement, there aie at least two known approaches of Kugo-Ojima (KO) ||T]] and 
Gribov-Zwanziger (GZ) [|2|]. In these approaches, the criteria for color confinement are attributed 
to the deep infrared behavior of specific Green functions. However, the Green functions depend on 
the gauge. Therefore, color confinement has been studied so far gauge by gauge, e.g., the Lorenz 
gauge. Coulomb gauge, the maximally Abelian (MA) gauge, etc. In the Lorenz-Landau gauge, 
especially, the KO criterion for color confinement is reduced to the infrared behavior of the ghost 
propagator 

Recent investigations show that gloun and ghost propagators in the most common Landau 
gauge are classified into two types according to their behaviors in the infrared region: 

• scaling solution (IR suppressed gluon propagator and enhanced ghost propagator) 
Schwinger-Dyson equation (SDE) [von Smekal, Hauck & Alkofer, 1997, 1998 [||]]... 

• decoupling solution (IR finite gluon propagator and tree-like ghost propagator) 
SDE [Boucaud et al.,2008 [|]] [Aguilar,Binosi, & Papavassihou,2008 [|]]... 

Lattice simulations [Bogolubsky et al., 2009 [|]] [Cucchieri & Mendes, 2008 0] [Sternbeck, 
von Smekal, Leinweber, Williams, 2007 ^]... 

The scaling solution fulfills the KO/GZ criterion, while this is not the case for the decoupling 
solution. The decoupling solution is supported by recent results of numerical simulations on the 
lattice with very large volumes [|6|, ^ |8|]. It is still under active debate to discriminate two different 
types of propagators. 

No one has found a gauge-independent criterion for color confinement! Therefore, quark 
confinement in the Wilsonian sense cannot be derived at present as a special case of these color 
confinement scenarios. In fact, the relationship between quark confinement and the Green functions 
in the infrared regime is not yet clarified. 

Nevertheless, it has been shown that both scaling and decoupling solutions exhibit quark con- 
finement and gluon confinement: 

• quark confinement, i.e. the vanishing of the Polyakov loop average at finite temperature T 
for < r < Tc below the critical temperature T^. 



functional renormalization group (FRG) [Marhauser & Pawlowski, 2008 QllO] [Braun, Gies 
and Pawlowski, 2010 [Kondo, 2010 [|l|]] 

gluon confinement, i.e., violation of reflection positivity was demonstrated 
for scaling solutions: SDE [Alkofer and von Smekal, 2001 [||]] 

for decoupling solutions: SDE,FRG [Fischer, Maas and Pawlowski, 2009 [^ ] Lattice simu- 
lations [Bowman et al., 2007] 
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The purpose of this talk is to discuss how quark confinement and gluon confinement are related to 



the infrared behavior of the gluon propagator in the Lorentz covariant gauge based on []13|]. 



(1) We derive a novel low-energy effective model of SU (2) Yang-Mills theory without fixing the 
original gauge symmetry. The resulting effective gluon propagator belongs to the Gribov- 
Stingl type, irrespective of the gauge choice: 

CQ + Cip^ +C2P^ 

(2) In MA gauge, we show that the model exhibits both quark confinement and gluon confine- 
ment simultaneously in the following sense: 

• quark confinement: The Wilson loop average satisfies the area law. 

• gluon confinement: A Schwinger function (Euclidean Green's functions) for the effective 
gluon propagator violates the reflection positivity. 

(3) However, for the effective gluon propagator to agree exactly with the Gribov-Stingl form 
Co 7^ 0, one must include either (a) a gauge-invariant nonlocal "mass term" or (b) a "mass 
term" that breaks nilpotency of the BRST symmetry. Otherwise, we have cq = 0. 

(4) We argue that quark and gluon confinement can be obtained even in the absence of such a 
mass term. 



In the follows, we consider only the SU (2) gauge group QlSj [19| , pO| ] and the extension to 



SU (N) based on [21] will be given elsewhere. 



2. Reformulating the Yang-Mills theory in terms of new variables 

(Step 1) We transform the original variables to the new variables n,c^, 

old variables : £^^{x) => new variables : (n'^ (x) , (x) , ^^^(x)), 
according to [T|, ||] 

n{x) =n'^{x)TA (A = 1,2,3) 
c^(x) = j2/^{x)-n{x), 

^^ix)=ig-'[D^[£/]nix),n{x)], (2.1) 

where n(x) is the Lie-algebra 5M(2)-valued field with a unit length, i.e., rt^{x)n^{x) = 1. The so- 
called color direction field n is obtained in advance as a functional of the original variable ^by 



solving the reduction condition [|T8||, e.g., [n(x),D^[^]Z)^[^]n(x)] = 0. 

If the original Yang -Mills field £/^{x) = £/^{x)Ta is decomposed into two pieces: =e^(x) = 
y^{x) + <f?r^(x), then the new variable = s^^ — .3^^ is a Lie-algebra i'M(2)-valued fields "^(x) = 
r^(x)rA (A = 1,2,3) given by 

y^ix) = Ci,{x)n{x) + ig-\n{x),d^n{x)]. (2.2) 

The variables ^(x) satisfy the properties: 
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(i) has the same gauge transformation as the original field =2/^, i.e., ^(x) — > Q.{x)Y^{x)D.{xy + 
ig'^Q.{x)d^Q.{xy and hence its field strength ^^v[f] ■= d^fy - dy^n - igi'fn,'^] trans- 
forms in the adjoint way: ^^v['^](-^) ^ Q.{x)^^v[y]a{xy , 

(ii) ^^iv['^] is proportional to n, i.e., ^^v['^] W := n(x)G^vW- 

Consequently, G^y = n- is gauge-invariant, since the field n is constructed so that it trans- 

forms as n(x) Q.{x)n{x)Q.{x)K Remarkably, G^y has the same form as the 't Hooft-Polyakov 
tensor for magnetic monopole: 

G^v = d^Cv - dyCi^ + ig^^n ■ [d^n,dvn]. (2.3) 

(Step 1') In order to obtain the dual effective theory for examining the dual superconductivity 
we introduce a gauge-invariant antisymmetric tensor field {*B)^y of rank 2 by inserting a 



(2.4) 



unity into the path-integral [22, E3, 24 1: 



1= J mtxp[- J d''x^{CB)^y-{an-^^y[r]-l5n-ig[^^,^y])y 



where * is the Hodge dual operation. Here (too many) parameters J, a, (5 aie introduced to see 
effects of each term. Putting 7 = is a simple way of reproducing the original theory without 
the antisymmetric tensor field B. When j3 = 7^^ = G and a = 0, indeed, (*B)^v is regarded as 
a collective field for the composite operator n ■ ig[3t^^, 3t^y\ with the propagator G obtainable in a 
self-consistent way [ ^ ] according to the Wilsonian renormalization group (RG) [^]. 

Then the Euclidean Yang-Mills Lagrangian ^ym[-2/] = li'^^vi-^])^ i^ rewritten and modified 

into 



^ym['^, ^ ,B\ 

= ^-^Gly + lCB)l, - f (*B),vG,v +^^^''Qfy^''' +^-^iig[^„^v])\ (2.5) 
where we have defined 



■.=S^''5^y + (2 + 7aj8)g£^^Cn^G^v - 7iSg£^^^«^(*S)^v, 
S^'':=-{Dp[r]Dp[y]f\ (2.6) 

with the covariant derivative in the adjoint representation with := y^Tc, (Tc)^^ = if^^^\ 
Df := d^S^^ - g/^^^rc = [dpl- igy.r. 



3. Deriving an effective model by eliminating high-energy modes 

(Step2') We identify with the "high-energy" mode in the range G [M^, A^] and proceed 
to integrate out the "high-energy" modes Here M is the infrared (IR) cutoff and A is the 
ultraviolet (UV) cutoff as the initial value for the Wilsonian RG. 

In the derivation of our effective model, we neglect quartic self-interactions among 3t^^, i.e., 
c^v])^- However, we can take into account an effect coming from the quartic interaction, 
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which influences our effective model. In fact, it is shown 1 2C , 35 1 that the quartic gluon interaction 
{ig[^^, ^v\Y among 3t^fi gluons can induce a contribution to the mass term 



1 2 



(3.1) 



through a gauge-invariant vacuum condensation of "mass dimension-2" (the BRST-invariant ver- 
sion was proposed in [|30|]), 



(3.2) 



which leads to the mass term ( jO) ) with ~ (^^^ (x) {x)) up to a numerical factor. This 
result is easily understood by a Hartree-Fock argument. This effect is included in the heat kernel 
calculation through the infrared regularization [|l^, 45]. 

The correlation functions for new variables have been computed on a lattice by numerical 



simulations using the Monte-Carlo method in [£^] based on [pq, [37[]. This justifies the identification 
of as the high-energy mode negligible in the low-energy regime below M ~ 1 .2GeV. Here the 
Landau gauge d^.s/^ = was adopted, since we need to fix the gauge to obtain the propagator or 
correlation functions. ' 

In these approximations, we can integrate out by the Gaussian integration and obtain a 
gauge-invariant low-energy effective action S'^^[f ,B\ 



^^^G2, + |(*B)2^-^(*B)p,Gp,] +ilndet(2^^-lndet5^« (3.3) 



where J = J d^x, the functional logarithmic determinant ^IndetQ^^ comes from integrating out 
the ^ field, and the last term comes from the FP-like determinant term [19] associated with the 
reduction condition [|I8|]. We can obtain (see a subsequent paper p5|] for details of calculations) 



llndetG^^,-lndet.- = /!^^ 

+ / ^^^(5A{(2 + 7«i8)Gp,-7i8(*S)p,})2 + 0(5VM^). (3.4) 



%-:^{{2 + wP)Gpa-rPCB)p^Y 



The gauge fixing is unnecessary in this calculation. Indeed, the resulting effective action ( |3.3| ) with 
(3.4) is manifestly gauge-invariant. This is one of main results. The correct RG j3 -function at the 
one-loop level j5{g) := = —big^ + 0{g^), b\ = ^/(47r)^ is reproduced in a gauge-invariant 

way when yajS = which follows from e.g. a = (a choice mentioned above) or 7 = (in the 
case of no B^v field). 

Thus we obtain the following effective action S'^^[G,B\ up to terms quadratic in the fields, 



G) + -{*B,[dQ + d,AYB) + {G, 



+ (3-5) 



'in the MA gauge, it has been shown in an analytical way that the off-diagonal gluon mass generation can 
follow from the off-diagonal gluon-ghost condensation of mass dimension 2, [g^ff),ff:= ^A"^A^" + aC"C" which has 
been proposed in [B0||. 
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where 

, 1 , ,,^2 g^ln^2-6(2 + y«^)^ 1 (2 + y«/3)^ 

/o-i + 7« ^^^^2 3 ' ■^i-(4;r)2M2 3 

We will see that the exact Gribov-Stingl form of the gluon propagator is obtained, if one 
introduces a gauge-invariant, but nonlocal "mass term": 

S'^^'[G] = ^{G,m^A-'G), (3.7) 
or, if one introduces a non gauge-invariant mass term 

S^^'[G] = ^{a,m^a). (3.8) 

for the gauge field a related to the field strength G by 

G = da, {5a = 0). (3.9) 

Even after taking specific gauges, the BRST invariance is also broken by including this mass term. 
However, we can modify the BRST such that the modified BRST is a symmetry of the Yang-Mills 
theory with the mass term at the cost of nilpotency. In other words, the requirement of nilpotency 
of the BRST excludes such a gluon mass term. 

4. The Gribov-Stingl form for gluon propagator 

We can integrate out B by the Gaussian integration. Then we obtain the effective action in 
terms of G: 

5ym[G] = ^(G,^ggG) (4.1) 
If we include the mass term (3.7), the inverse effective propagators for the field strength G reads 



i^G^ =[m^A-'+fo + fiA] - [do + diA]-\ho + hAf 

Jdo + diA][m^A-'+fo + fiA]-[ho + hA]^ ^4 2) 
do + diA 

Then we obtain the effective propagator for the field a defined by G = da: 

^GG-^ ^aa, ^aa " ^^^^^^ , (4.3) 

where 

Co =m^dQ, ci = Jo/o - + ^di , C2 = dofi + /o<ii - 2/jo/ii . (4.4) 
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We observe that the effective gluon propagator 

%a{p) 



^aa has the Gribov-Stingl form when cq 7^ 0: 



(4.5) 



Thus, we have found that the effective propagator has the Gribov-Stingl form []14[]. Note that di 
comes from the induced kinetic term for the B field which was introduced in the beginning as an 
auxiliary field without the kinetic term. 



5. Converting the Wilson loop to the surface-integral 



(Step 3) We use a non-Abelian Stokes theorem [ ]25| , [26| , [27| ] to rewrite a non-Abehan Wilson 
loop operator 



Wc[s^] :=tr 



^ exp i ig j dx^s^^ (x) 



Ar(l), 



into the area-integral over the surface £ (5£ = C): 

1 



WcK]= j diizi^)exp 



(5.1) 



(5.2) 



where d\l is an invariant measure on SU (2) normalized as / d\i(^^^ = 1, S SU (2). In the two- 
form G := ^G^v{x)dx^ Adx^, G^y agrees with the field strength ( |2.3| ) under the identification of 
the color field n{x) with a normalized traceless field (See also [^]) 



(5.3) 



Using the vorticity tensor &^ with the support on the surface £ whose boundary is the loop C: 

d^S^\x{c>))d'^{x-x{a)), (5.4) 



the surface integral is cast into the volume integral and the Wilson loop operator is rewritten as 



(0z,G) = I J^x^Qf (x)G^vW, (5.5) 



where (•, •) is the inner product for two differential forms. 



6. Calculating the Wilson loop average to show area law: quark confinement 

(Step 4) We proceed to evaluate the Wilson loop average W(C) = (Wc[^])ym by using the 
effective action Sfj^[G,B], i.e., {Wc[^])ym ^ (Wc[^/])ym with the aid of (|^). 

We demonstrate that the simplest way to obtain the area law is to use the low-energy effective 
action 5^[G,B] retained up to terms quadratic and bilinear in G and B. 

In what follows, we take the unitary-like gauge 

n^{x) = 5a3, (6.1) 
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which reproduces the same effect as taking the MA gauge [^] in the original Yang-Mills theory. 
In this gauge, (x) reduces to the off-diagonal component (x) (a = 1,2), while (x) reduces 
to the diagonal one A^(a;) = a^{x), i.e., ^^(x) = £/^{x)5Aa, ^C-^) = ■^^{x)5a3 = c^(x)5a3. The 
gauge ( |6.1[ ) forces the color field at each spacetime point to take the same direction by gauge 



rotations. Hence the field G given by ([2.3|), 
G 



dyC^ , 



H^y ■=ig ^n-[d^n,dyn] 



(6.2) 



contains singularities (of hedgehog type) similar to the Dirac magnetic monopole after taking the 
gauge (^J]). If we do not fix the gauge, such a contribution is contained also in the part ig^^n ■ 



[d^n,dyn] to make a gauge-invariant combination G^y, see [ pq , 37]. Consequently, the Bianchi 
identity for G is violated, 



5*G = *dG = **ddc + *dH = 5*H ^0, 



(6.3) 



even if ddc = 0. Here d denotes the exterior differential and 5 the codifferential. There is no 
well-defined one-form h such that H = dh. Thus we obtain a nontrivial gauge-invariant magnetic 
monopole current defined by 

k:=5*G. (6.4) 

By integrating out the B field, we obtain the effective action 5^[G]. Then we find that the 
effective propagator has the Gribov-Stingl form: 



CQ + Cip^ +C2P^ 



(6.5) 



where cq ■ 



■ rrp-, ci 



1 + 



3 (4;r)2M2^ 



'^2=(^ps|i[(2 + 7«i8)2 + (l + 7«')7i8' + 2(2 + 7«i8)7«i8]/3, 



and d 



3 {47Z) 



^-gy. The precise values of the parameters m,Y,a,p and M are to be determined 



by the functional RG [29] following [12], which is a subject of future study. 



In the unitary-like gauge (|6JJ) the Wilson loop operator is reduced to 

Wc[F]=exp 

Then the Wilson loop average W(C) is evaluated by integrating out G = da: 

1 











= exp 




2 .Jl,:dT,=C 





W(C) =exp 



(6.6) 



(6.7) 



where &gg = 4^aa and its Fourier transform obeys S^Goip) = P^^aaip)- 

For concreteness, we choose @y. for a planar surface bounded by a rectangular loop C with 
side lengths T and R in the X3 — X4 plane. Then we find that the Wilson loop average has the area 
law for large R and T: 

W(C) ~exp[-a/?r], (6.8) 
with the string tension given by the formula: 

^ 8^ Jp\+pl<M'- ' (271) 



^^!^GGiPuP2,0,0)>0, 



(6.9) 
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Figure 1: (Left panel) The propagator of the diagonal gluon in MA gauge as a function of the momentum 
P — \f^- (Right panel) The dressing function of the diagonal gluon in MA gauge, or the integrand of the 

\p\ „ fiinr-tirvn r\f n — a / n2 



formula for the string tension, — , i^"^' , as a function of p = v n 

where the momentum integration is restricted to the two-dimensional momentum space (the di- 
mensional reduction by two pi]]) and is cutoff at M which is the upper limit of the low-energy 
effective model being meaningful. A positive and finite string tension < a < oo follows from 
the condition of no real poles in the effective gluon propagator ^gg(/') in the Euclidean region, 
< ^gg{p) = P^^aa{p) < °°, which is connected to the gluon confinement shown below. This is 
another main result. 



According to numerical simulations in MA gauge |^, ^ ^J], the diagonal gluon propagator 
is well fitted to the form dJ), see Fig. 1: e.g. @| gives cq = 0.064(2)GeV2, ci = 0.125(9), 



C2 = 0.197 (9)GeV-2, di = 0.13(l)GeV-2, andM ~ 0.97GeV, where M is the mass of off-diagonal 
gluons obtained in the MA gauge. This value of M is a little bit smaller than the values of other 



groups [ |42| , |43| |. This indeed leads to a good estimate for the string tension 

a~(0.5GeV)2, (6.10) 



according to ( |6!9| ) for a(/x) = g^(/i)/(47r) ~ 1.0 at /i = M. The next task is to study how the 
results are sensitive to the deep infrared behavior of the diagonal gluon propagator ( |63| ) and the 
actual value of M for the off-diagonal gluon propagator. 

The Gribov-Stingl form is obtained only when cq^O (i.e., m^O) and di y^O (B^y is included). 
Even in the limit — )• (co — )• 0), the area law can survive according to (^^, provided that 
^gg(p) remains positive and finite: ^ooip) ci ' ^^^^^ ^aa(p) behaves unexpectedly as 

^aaip) p2(c^^'c^p2-) - Hence, we argue that it does not matter to quark confinement whether m = 
or m 7^ 0. 

7. Calculating the Schwinger function to show positivity violation: gluon 
confinement 

(Step5) The positivity violation is examined. We consider the Schwinger function defined by 

A{t) := J d'xe-'P-^&{t,x)\p=o = f^^'"'"^^P = ^^P'^)- ^^.l) 
The Euclidean propagator ^{p) in momentum space has a spectral representation, 

^{p) = HdK^ > ^ A(0 = rdKp{K^)e-^' > 0. (7.2) 

Jq + Jo 
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Figure 2: The Schwinger function A{t) calculated from the diagonal gluon propagator of the Gribov-Stingl 
type obtained in MA gauge (left panel) < f < 5, (middle panel) < f < 12, (right panel) 5 < f < 15. 

If A(f) is found to be non-positive, p(k:^) cannot be a positive spectral function. The corresponding 
states cannot appear in the physical particle spectrum: they are confined. 
1) For the free massive propagator, A{t) is positive for any t, 

Therefore, we find no positivity violation as expected. This case corresponds to p(K"^) = 5(k"^ — 



m 



i^5(jc-m)>0. 



2m 



2) We consider the propagator of the Gribov-Stingl type in Euclidean space, 

^{p) = —^^^^4^ — I, P^>0, co,cuC2,do,di£R. (7.4) 

In the case of C2 = 0, there is no positivity violation, as far as co/ci > 0. In the case of C2 / 0, 
^aa{p) has a pair of complex conjugate poles at =z and p^ = z*, z :=x+iy,x:= —ci/{2c2),y-= 

co/c2 - (ci/(2c2))^. We find that the Schwinger function A(0 := /+r ^e'P*'^^{p = 0,p^) is 
oscillatory in t and is negative over finite intervals in the EucHdean time t > (See Fig. 2): 

AU) = — i — - — -f'"'W'^'^'"'''[cos(f|z|i/2cos(p-(p) + Ji|z|cos(f|z|i/2cos(p + (p)l, (7.5) 

2c2|z|-*/^sin(2<p) 



wherez= with |z| = (co/c2)'^^, cos (29) = -^cy{4coC2), and sin (2 9) = ^J 1 -c^/(4coC2) 

Therefore, the reflection positivity is violated for the gluon propagator (|6.5[), as long as 



0<-^<l, (7.6) 

4coC2 

irrespective of di . When co = (or m = 0), 

Hence, the special case cq = also violates the positivity, if ci > and C2 > 0. Thus the diagonal 
gluon in the MA gauge can be confined. 
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8. Summary 

In summary, we have discussed how to obtain a low-energy effective model of the SU (2) 
Yang-Mills theory without fixing the original gauge symmetry. It is remarkable that the effective 
model respects the SU (2) gauge invariance of the original Yang-Mills theory, which allows one to 
take any gauge fixing in computing physical quantities of interest in the low-energy region. The 
resulting effective gluon propagator belongs to the Gribov-Stingl type, irrespective of the gauge 
choice This is a universal aspect obtained independently of the choice of gauge fixing condition. 

In MA gauge, especially, we have demonstrated that the model exhibits both quark confine- 
ment and gluon confinement simultaneously in the sense that the Wilson loop average satisfies the 
area law (i.e., the linear quark-antiquark potential) and that the Schwinger function violates re- 
flection positivity. Moreover, we have given a formula for the string tension based on the gluon 
propagator of the gauge-invariant field strength G^y It gives a good estimate for the string tension. 

However, for the effective gluon propagator to agree exactly with the Gribov-Stingl form, we 
need to introduce (i) a gauge-invariant, but nonlocal mass term or (ii) a mass term that breaks 
nilpotency of the BRST symmetry. Otherwise, we have cq = 0. We argued that both quark and 
gluon confinement can be obtained even in the absence of such a mass term, cq = 0. More results 
and full details will be given in a subsequent paper [p^]. 
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